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Low-momentum nucleon-nucleon interactions are derived within the framework of a unitary- 
transformation theory, starting with realistic nucleon-nucleon interactions. A cutoff momentum 
A is introduced to specify a border between the low- and high-momentum spaces. By the Faddeev- 
Yakubovsky calculations the low-momentum interactions are investigated with respect to the de- 
pendence of ground-state energies of 3 H and 4 He on the parameter A. It is found that we need the 
momentum cutoff parameter A > 5 frn -1 in order to reproduce satisfactorily the exact values of the 
binding energies for 3 H and 4 He. The calculation with A = 2 fm _1 recommended by Bogner et al. 
leads to considerable overbinding at least for the few-nucleon systems. 

PACS numbers: 21.30.-x, 21.45.+V, 27.10.+h 



I. INTRODUCTION 

One of the fundamental problems in nuclear structure 
calculations is to describe nuclear properties, starting 
with realistic nucleon-nucleon (NN) interactions. How- 
ever, since this kind of interaction has a repulsive core at 
a short distance, one has to derive an effective interaction 
in a model space from the realistic interaction, except for 
the case of precise few-nucleon structure calculations. 

Recently, Bogner et al. have proposed a low- 
momentum nucleon-nucleon (LMNN) interaction which 
is constructed in momentum space for the two-nucleon 
system from a realistic nucleon-nucleon interaction, using 
conventional effective interaction techniques or renormal- 
ization group ones [| . In the construction of the LMNN 
interaction a cutoff momentum A is introduced to specify 
a border between the low- and high-momentum spaces. 

The LMNN interaction is constructed in order to ac- 
count for the short-range correlations of the two nucleons 
interacting in the vacuum. So the question is to what 
extent the obtained LMNN interaction is a good approx- 
imation also for describing correlation of nucleons inter- 
acting in a nuclear many-body medium. The medium 
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effect could appear through the single-particle potential 
and three-or-more-body correlations. From a practical 
point of view it is of high interest to explore the sensitiv- 
ity of calculated results to the cutoff momentum A. 

Bogner et al. have constructed their LMNN interac- 
tion in a way to conserve in the low-momentum region 
not only the on-shell properties of the original interaction 
(i.e. phase shifts and the deuteron binding energy) but 
also the half-on-shell T matrix pj. They found that the 
LMNN interactions for A = 2.1 fm _1 corresponding to 
Slab ~ 350 MeV become nearly universal, not (or only 
weakly) dependent on the choice of realistic interactions 
employed. They suggested to use their LMNN interac- 
tion directly in nuclear structure calculations, such as 
the shell-model 0] and the Hartree-Fock calculations . 
They claimed that for the momentum cutoff in the vicin- 
ity of A = 2.1 fm _1 the calculated low- lying spectra for 
18 0, 1 Te, and 135 I are in good agreement with the ex- 
perimental data and depend weakly on A. Their results 
have been found to agree with the data as good or even 
slightly better than the results based upon the G matrix 
which is constructed by taking into account the short- 
range correlations, the Pauli blocking effect, and the state 
dependence for each nucleus. 

Kuckei et al. 4] investigated the nuclear matter and 
the closed-shell nucleus 16 O by using the LMNN inter- 
action obtained by Bogner et al. They concluded that 
the LMNN interaction can be a very useful tool for low 
energy nuclear structure calculations, and one should be 
cautious if the observable of interest is sensitive to the 
single-particle spectrum at energies above the cutoff mo- 
mentum. 
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However, one had better compare with the exact so- 
lutions once for all. In cases of three- and four-nucleon 
systems we can directly make such a comparison by solv- 
ing the Faddeev-Yakubovsky equations 



We investigate the LMNN interaction by means of uni- 
tary transformation using two independent (but equiva- 
lent) approaches. One approach is based on the uni- 
tary transformation of the Okubo form 6] in which a 
LMNN interaction is obtained from the scattering am- 
plitude in momentum space 0, H, 0- Another is the 
unitary-model-operator approach (UMOA) 0, 
contrast to the G-matrix theory the UMOA leads to 
an energy-independent and Hermitian effective interac- 
tion in a many-body system. Contrary to the LMNN of 
Ref. pi, we will not require conservation of the half-on- 
shell T matrix, which does not represent an observable 
quantity. Only low-momentum NN observables such as 
the on-shell T matrix, phase shifts and binding energies 
are guaranteed to remain unchanged under an unitary 
transformation. In principle, one could achieve the equiv- 
alence of the half-on-shell T matrix by performing an 
additional unitary transformation in the low-momentum 
space. We, however, refrain from doing that since we do 
not see any conceptual or practical advantage in requir- 
ing the equivalence of the half-on-shell T matrix. 



In the present study we first apply the above- 
mentioned methods to the two-nucleon system in momen- 
tum space to construct the LMNN interaction. Although 
both approaches are based on the same idea of an uni- 
tary transformation, the calculation procedures for de- 
riving the LMNN interaction are independent and quite 
different from each other. We calculate selected prop- 
erties of the two-nucleon system using both schemes to 
confirm the numerical accuracy. This cross check is use- 
ful to ensure the reliability of the numerical calculations. 
Secondly, we investigate the A dependence in structure 
calculations of few-nucleon systems, where (numerically) 
exact calculations can be performed 0, , and discuss 
the validity of the LMNN interaction by comparing the 
obtained results with the exact values. 



This paper is organized as follows. In Sec. [Q] af- 
ter the basic formulation of the unitary transformation 
is presented, the two methods are given with empha- 
sis on the different calculation procedures for deriving 
the LMNN interaction. In Sec. EH LMNN interac- 
tions are constructed using both methods from realistic 
nucleon-nucleon interactions such as the CD-Bonn prl ] 
and the Nijm-I 15] potentials. Then, the Faddecv and 
Yakubovsky equations are solved using the LMNN in- 
teractions for various values of the cutoff parameter A. 
Finally, we summarize our results in Sec. IIVI 



II. UNITARY TRANSFORMATION OF THE 
HAMILTONIAN FOR THE TWO-NUCLEON 
SYSTEM IN MOMENTUM SPACE 

We consider a quantum mechanical system described 
by a Hamiltonian H . The Schrodinger equation reads 



P* = EV. 



(1) 



Introducing a unitary transformation U with UIP = 1 
we obtain a transformed Schrodinger equation 



H'V = EM' 



(2) 



with the transformed Hamiltonian and state, H' = 
U^HU and = U'^f, respectively. We introduce the 
concept of the effective Hamiltonian by means of the par- 
tition technique. The Hilbert space is divided into two 
subspaces, the model space (P space) and its complement 
(Q space), such that the Schrodinger equation becomes 
a 2 x 2 block matrix equation 



PH'P PH'Q 
QH'P QH'Q 





' p<s>>~ 


= E 


' P*' ' 









(3) 



Here P and Q are the projection operators of a state 
onto the model space and its complement, respectively, 
and they satisfy P + Q = 1, P 2 = P, Q 2 = Q, and 
PQ = QP = 0. The Q-space state is easily eliminated to 
produce the projected Schrodinger equation 

[ PH'P + PH'Q^j^QH'P ] (P*') = E{PV). (4) 

In general, the effective Hamiltonian which is given in 
parentheses depends on the energy E to be determined. 
However, if the decoupling equation 



QH'P = 



(5) 



is satisfied, then we have the equation for the energy 
independent effective Hamiltonian in the P space 



PH'P(Pm') = E(P^'). 
A unitary transformation can be parametrized as 



(6) 



U 



P(l+«t W )-V2p -Pcjt(l + c ^t)-l/2( 

Qw(i + (Aj)- 1/2 p Q(l + ^t)-l/2Q 



(7) 

and the wave operator o> satisfies the condition u> — QluP. 
Equation J7J is well-known as the Okubo form . Notice 
that the unitary transformation given in Eq. JJJ is by no 
means unique: in fact one can construct infinitely many 
different unitary transformations which decouple the P 
and Q subspaces. For example, performing subsequently 
any additional transformation, which is unitary in the 
P subspace, one would get a different LMNN interac- 
tion . The transformation in Eq. (7J depends only on 
the operator to which mixes the P and Q subspaces and 
is in some sense "the minimal possible" unitary trans- 
formation. For more discussion the reader is referred to 
Ref. @. 
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In the present study, the above unitary transformation 
is used in two different methods to derive LMNN inter- 
actions. In the following sections, we shall give details on 
the two methods. 



A. Method-1 

Consider a momentum-space Hamiltonian for the two- 
nucleon system of the form 



H(p,p') = H (p) 8{p~ p') + V{p,p') , 



(8) 



where Ho(p) — p 2 /(2M) with the reduced mass M 
stands for the kinetic energy, and V(p,p') is the bare 
two-body interaction. Our aim is to decouple the low- 
and high- momentum components of this two-nucleon po- 
tential using the method of unitary transformation. To 
achieve that, we introduce the projection operators 



P = J d 3 P \p)(p \ , \p\ < A , 



Q 



d 3 q\q)(q\, \q \ > A 



(9) 



where A is a momentum cutoff whose value will be spec- 
ified later, and P (Q) is a projection operator onto 
low- (high-) momentum states. Based on the unitary- 
transformation operator given in Eq. 0, the effect lve 
Hamiltonian in the P space takes the form 

PH'P = P(l + lu^uj)^ 1 / 2 (H + uj^H + Hlu + lu^Hu) 
x {1+^lu)-^ 2 P. (10) 

This interaction is by its very construction Hermitian [(| . 

The requirement of decoupling the two spaces leads to 
the following nonlinear integral equation for the operator 

V(q,p) - J dYuWW,?) 

+ / d 3 q'V(q,q')u;(q',p) 



d 3 p'd 3 q'ij(q,p')V(p',q')ij(q',p) 

= (E p -E q )u(q,p), (11) 

where we have denoted by p (q) a momentum of the P 
space (Q space) state and by E p (E q ) the kinetic energy 
E p =p 2 /(2M) (E q =q 2 /(2M)). 

Alternatively, one can determine the operator oj from 
the following linear equation, as exhibited in Refs. 0, @ : 



of the Lippmann-Schwinger (LS) type since the position 
of the pole E p in the integration over p' is not fixed but 
moves with p. In solving the integral equation Eq. 1)12(1 , 
the second argument p in u> varies, whereas the first one 
q is a parameter. 

In this study we have used the linearised equation 
Eq. (|12[1 to project out high-momentum components from 
the realistic potentials. Since the projection operators P 
and Q do not carry any angular dependence, the integral 
equation Eq. H12[l can be solved for each partial wave 
independently. In the partial wave decomposed form it 
reads 



Tu!(q,p,E q ) 



E n — E„ 



(13) 



E 



Jp 



, u s 2{q,p>)T%(jp>,p,E p ) 
p z dp ' 



E p — E- 



p' 



where V$(q,p) = (lsj,q\V\l'sj,p) and u)${q,p) = 
(Isj, q\u>\l'sj,p). In the uncoupled case I is conserved and 
equals j. In the coupled cases it takes the values I = j± 1. 

Equation (|12|) and, consequently, also Eq. (|13|) have a 
so-called moving singularity, which makes it more diffi- 
cult to handle than the LS equation. Indeed, one has to 
discretize both p and p' points in Eq. (fT5f> . This does 
not necessarily allow to solve Eq. I|13l) . The difference 
E p — Ep' can be exactly zero since p and p' belong now 
to the same set of quadrature points. Thus, one cannot 
calculate the principal value integral in the same manner 
as for the LS equation. To solve this equation we have 
used a method proposed by Glockle et al. 9]. 

Another problem arising in solving the linear equa- 
tion is caused by the fact that the driving term 
Tip (q, p, E p ) / (E p — E q ) becomes very large when p = \p \ 
and q = \q \ go to A. Consequently, the equation be- 
comes ill-defined. To handle this problem we have reg- 
ularized this equation by multiplying the original poten- 
tial V(k', k) with some smooth functions f(k') and f(k) 
which are zero in a narrow neighborhood of the points 
k' = A and k = A ^tJ- The precise form of this regular- 
ization does, in fact, not matter 

Having determined the operator u>, one can calculate 
the effective Hamiltonian in the P space according to 
Eq. (TD}. To evaluate the operator P(l + uj^uj)^ 1 / 2 en- 
tering this equation we first diagonalize the operator 
(1+o^cj) and then take the square root. Finally, the effec- 
tive potential can be found by subtracting the (not trans- 
formed) kinetic energy term from the effective Hamilto- 
nian. 



= T{q 7 p,E p ) _ f bj(q,p') T(p',p,E p ) 
E p — E q J E p — E p ' + ie 



Method-2 



Here the integration over p' goes from to A. Conse- In Eq. J7J, the unitary-transformation operator U has 
quently, the dynamical input in this method is the T ma- been given by the block form with respect to the projec- 
trix T(pi,p2, E P2 ). Note that this is not a usual equation tion operators P and Q. We notice here that the operator 
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U can also be written in a more compact form as ^| 

U ={1 + uj-uj^){1+ujuj^ +uj ] uj)- 1/2 . (14) 

Using the above operator U, in general, the effective in- 
teraction V for a many-miclcon system is defined through 



V = U- 1 (H + V)U - Ho 



(15) 



where Hq is the kinetic energy of the constituent nuclc- 
ons in the nuclear system, and V is the bare two-body 
interaction between the nucleons. Here we apply thus 
defined effective interaction to the two-nucleon problem. 
In that case Hq becomes the relative kinetic energy of 
the two nucleons, and V the bare two-body interaction 
between the two nucleons. Then the LMNN interaction 
Mow-fc of interest in the present work is given by 



Via 



= P 



low 



-kVPl 



(16) 



where -Pi ow -fe is the projection operator onto the low- 
momentum space for relative two-body states, and is the 
same as P in Eq. @ of the method- 1. In order to ob- 
tain Mow-fc in the form of the matrix elements using the 
plane-wave basis states, we shall present in the following 
a procedure for the numerical solution. 

We first consider an eigenvalue equation for the relative 
motion of a two-nucleon system as 



(H + V)\* r , 



(17) 



The above equation is written also in an integral form 
concerning relative momenta k and k! as 



(k'\H + V\k)(k\^ n )k 2 dk = E n {k'\V n ), (18) 



where 



\k)(k\k 2 dk = 1. 



(19) 



We here make an approximation for Eq. (|18|) in a nu- 
merical integral form by introducing the adequate inte- 
gral mesh points ki and kj and discretizing k — > kj and 

J2(i*i\Ho + V\k 3 ){k^ n ) = E n (ki\* n ), (20) 



where \~k~i) and \kj) represent the plane- wave basis states 
in the matrices. Those grids are characterized by the 
mesh points ki and kj. Thus, is defined as 



\h) = hVwl\h) 



(21) 



with the plane- wave states \ki) and the weight factors for 
the numerical integral \JWi- They are normalized as 



(ki | kj 



(22) 



We note here that the eigenvectors \^ n ) in Eq. (|17fl can 
be expressed as |^ n ) = \(j) n )+uj\(j) n ) in terms of the opera- 
tor lj and the P-space components \(f> n ) = P\^ n ). Thus, 
the formal solution of ui is given by w = J2 n \^n)(<t>n\ 
with the biorthogonal state (<j> n \ of | <j> n ). In order to ob- 
tain the matrix elements of u, we first solve Eq. (|20|) by 
diagonalizing the matrix elements, using the basis states 
\ki). Then, the matrix elements to for the basis states 
\kp) in the P space and \k q ) in the Q space are obtained 
as 



(A,M*p) = Y,(k q \Q\* n )(4> n \P\k P ), 



(23) 



where d is the number of the basis states (the inte- 
gral points) in the P space. As shown in Eq. © 
the P space (low-momentum space) and Q space (high- 
momentum space) are defined with a cutoff momentum 
A as < k < A and A < k < oo, respectively. The 
bra states (4> n \ are obtained through the matrix inver- 
sion as [(<^n|fcp)] = [(k p > l^n)]" 1 and satisfy the relations 



Efc TJ (^|fcp)(fcp|0„') 



5k p ,.k p - It should be noted that the solution w given in 
Eq. 12: il) is ambiguous in the sense that how to choose the 
set of eigenvectors {\^ n ),n — 1,2, ...,d} is not unique. 
We here select {l^n)} so that they have the largest P- 

space overlaps O n — Yl%=i \(^n\P\ki)\ 2 among all the 
eigenstates in Eq. (|20|) . As we will show later, the nu- 
merical calculation shows that this selection of {l^n)} 
leads to the same solution obtained from method- 1. 

In order to obtain the LMNN interaction (P-space ef- 
fective interaction) , we introduce the eigenvalue equation 
for ui'uj in the P space as 



(24) 



Using the solutions to the above equation, the LMNN 
interaction of a Hermitian type is given by |18l [l9j 

0'a\Vlow-k\i>f3) 



y/(l+»lM a \R\ll>p} + ^(l+M^al^l^) 

v / 0+74) + ^/(i + /4) 



where 



R = P(V + Vuj)P 



(25) 



(26) 



is a low-momentum (effective) interaction of a non- 
Hermitian type. Finally, the matrix elements of the 
LMNN interaction using the plane- wave basis states \ki) 
and \ks) are obtained as 



{ki\V\ oyi -k\kj) = 



Ea,/i I V'a) (tp a \V\ ow - k \tpp) (lpp\kj) 

k~k~JWW 3 



( 2? ) 

By an interpolation technique the elements of the poten- 
tial are prepared at arbitrary momenta in the P space. 
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FIG. 1: Comparison of LMNN interactions from the CD- 
Bonn potential in the case of A = 2.0 fin -1 . The diagonal 
matrix elements (k\ Vi ow _fe|fc) for the 1 So (a) and 3 5i- 3 Di (b) 
partial waves are shown. The lines depict the method-1 (solid) 
and the method-2 (long-dashed) results. Because they coin- 
cide very well, one cannot distinguish both lines by the eye. 



III. RESULTS AND DISCUSSION 

As mentioned in the previous section, we have two 
different methods based on the unitary transformation. 
Here we make the cross check using both methods. 

In Fig. n the diagonal matrix elements of the LMNN 
interactions for the neutron-proton 1 Sq and 3 Si- 3 Di 
channels using the CD-Bonn potential l3| are shown 
in the case of A = 2.0 fm~ 1 . In order to see the off- 
diagonal matrix elements, we also illustrate (fc|14ow-fe|2fc) 
in Fig. [21 One can see that the results obtained by the 
two methods are almost the same for both the diagonal 
and non-diagonal matrix elements within 3-4 digits using 
typically 100 integral grid points. Having obtained the 
same results with high precision using two very different 
methods, we have confidence in our numerical results. 

In Fig. |3| we show the same matrix elements of the 
original CD-Bonn potential and the LMNN interaction 
for the sake of comparison. The LMNN potential is 
very different from the original one. Nevertheless, the 
scattering phase shifts and the mixing parameter below 
E lah = 300 MeV for the and 3 S'i- 3 £)i channels, as 



FIG. 2: Comparison of LMNN interactions for the off- 
diagonal elements (k\Vi ow -k\2k) . Description is the same as 
in Fig. □ 



shown in Fig. reproduce exactly the ones obtained 
from the original interaction. This has been also shown 
by Bogner et dl. 

We now regard deuteron properties. In the light of the 
unitary transformation for the two-nucleon system, all 
the calculated deuteron binding energies for various val- 
ues of A must reproduce the exact value using the original 
interaction. In TableQJ calculated deuteron binding ener- 
gies for various A using the CD-Bonn [14| and the Nijm- 
I potentials are tabulated together with deuteron in- 
state probabilities. Note that the -D-state probability is 
not observable [20|. The values in the last row are the 
exact values quoted from the original papers of the CD- 
Bonn and the Nijm-I potentials. Indeed, we see that the 
calculated binding energy for each A reproduces the ex- 
act value with high accuracy. However, as for the -D-state 
probability, the difference between the results using the 
LMNN interaction and the exact value becomes larger 
as the value of A becomes smaller. If one uses the cor- 
responding effective operator in calculating the D-state 
probability, i.e. the unitarily transformed projection op- 
erator onto the D state, then one would reproduce the 
original value for this quantity. 

As shown in Table [I] for the deuteron, the difference 
of the wave functions does not affect the binding energy 
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FIG. 3: Comparison of the LMNN interaction in the case of 
A = 2.0 far 1 and the original CD-Bonn potential for the So 
(a) and 3 S\- 3 Di (b) partial waves. The solid and dashed lines 
depict the diagonal matrix elements of the LMNN interaction 
and the original CD-Bonn potential, respectively. 
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channels below i?i a b = 300 MeV. Because the lines from the 
LMNN interactions and the original CD-Bonn potential coin- 
cide, one cannot distinguish both lines by the eye. In Fig. (b) 
the upper, middle, and lower lines depict the 3 5i phase shift, 
the mixing parameter ei, and the 3 Di phase shift, respec- 
tively. 



TABLE I: Calculated binding energies Z?b(MeV) and D-state 
probabilities Pd(%) of the deuteron for various values of A. 
The values in the last row are those quoted from the original 
papers of the CD-Bonn and the Nijm-I potentials. 



CD Bonn 



Nijm I 



A(fm~ 1 ) 


Sb(MeV) 


Pd(%) 


£ b (MeV) 


Pd(%) 


1.0 


-2.224576 


1.21 


-2.224575 


1.24 


2.0 


-2.224576 


3.55 


-2.224575 


3.83 


3.0 


-2.224576 


4.55 


-2.224575 


5.12 


4.0 


-2.224576 


4.79 


-2.224575 


5.53 


5.0 


-2.224576 


4.83 


-2.224575 


5.64 


6.0 


-2.224576 


4.83 


-2.224575 


5.66 


7.0 


-2.224576 


4.83 


-2.224575 


5.66 


quoted 


-2.224575 


4.83 


-2.224575 


5.66 



since we perform the unitary transformation for the two- 
nucleon system. However, if we apply the LMNN in- 
teraction to the calculation of the ground-state energies 
of many-body systems, the situation will change. This 
is because the unitary transformation in the two-nucleon 
Hilbert space is not unitary any more in the Hilbert space 
of three and more nucleons. As a consequence, calcu- 



lated binding energies will depend on A. In order to 
examine the A dependence, we perform the Faddeev and 
Yakubovsky calculations for few-nucleon systems. Re- 
cent precise calculations for few-nucleon systems are re- 
viewed in Refs. [3 [3. 

Figure EJa) exhibits the energy shift A£"b = Eb(A) — 
£"b(oo) from the ground-state energy £b(oo) of 3 H as a 
function of A based on the CD-Bonn potential (solid) 
and the Nijm-I one (short-dashed) by a 34-channel Fad- 
deev calculation. In the present study, only the neutron- 
proton interaction is used for all the channels for sim- 
plicity. The exact value -Eb(oo) using the original po- 
tential CD Bonn (Nijm I) is -8.25 (-8.01) MeV. The 
long-dashed line depicts the Faddeev calculation using 
the original interaction, where the high-momentum com- 
ponents beyond A are simply truncated. The numerical 
stability is lost within the area A < 1 fin" 1 . For the 
case of the CD-Bonn potential, we need A > 8 fin -1 
to reach the exact value if the accuracy of 100 keV is 
required for the case of simple truncation. This situ- 
ation is largely improved if we use the LMNN interac- 
tion. Even if we require the accuracy of 1 keV, we do 
not need the high-momentum components beyond A ~ 8 
fin -1 . However, it should be noted that the results us- 
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FIG. 5: Energy shifts AEb of 3 H (a) and 4 He (b) as a function 
of A using the LMNN interaction (solid) from the CD-Bonn 
potential. The long-dashed lines in (a) and (b) are plotted for 
the case of simple momentum cutoff calculations. The short- 
dashed line in (a) depicts the result for the LMNN interaction 
from the Nijm-I potential. 



ing the LMNN interaction for the values smaller than 
A ~ 5 fm -1 vary considerably, and there occurs the en- 
ergy minimum around A ~ 2 fin . The cutoff value, 
which produces the minimal value of A£"b , is close to the 
value proposed by Bogner et al. [[J. 

We note here that the wave function using the LMNN 
interactions is very close to the true wave function if the 
Jacobi momentum set (p, q) of the three-nucleon system 
is smaller than the adopted value of A when we take 
A > 4 fin- 1 . 

By solving the Yakubovsky equations we can calculate 
the binding energy of 4 He. A similar tendency of the A 
dependence for 3 H can be also seen in the results for 4 He. 
In Fig.[SJb), the A dependence of the energy shift AEt, 
of the ground state of 4 He is illustrated. We demonstrate 



only the case of the CD-Bonn potential. Here we adopt 
the S-wave (5 + 5-channel) approximation and do not in- 
clude the Coulomb force for simplicity We also use 
only the neutron-proton interaction for all the channels. 
The exact value E\>(oo) using the original CD-Bonn po- 
tential in the above-mentioned approximations is —27.74 
MeV. The shape of the curve is similar to that for 3 H up 
to A ~ 2 fm -1 . The numerical instability already starts 
at A < 2 fin -1 . Here it seems to conserve the strong cor- 
relation relation between AE h ( 4 Ue) and AS b ( 3 H) 
where the ratio is about 5 in the region A > 3 fm -1 . 

In the case of 3 H, the energy shift AE\> is about 750 
keV (1 MeV in case of Nijm-I potential) at the minimum 
point. On the other hand, the difference amounts to 
about 3 MeV for 4 He at A = 2.0 fin -1 . We remark 
that the minimal value of the cutoff A which leads only 
to small deviations of the 3 H and 4 He binding energies 
from their exact values is A ~ 5 fin -1 . As far as the 
excitation spectra of low-lying states from the ground 
state are concerned, the A dependence at A ~ 2 fin -1 
may be weak and the LMNN interactions could be useful 
as has been shown in the shell- model calculations 0. 

It is noted that for LMNN interactions considered as 
a purely computational tool, which allows to get rid of 
the hard core and thus enables many-body calculations 
based on realistic NN potentials without using the G- 
matrix formalism, we recommend the value of the cutoff 
A of the order of (or higher than) ~ 5 fin -1 . One then 
ensures that the binding energies do not deviate signifi- 
cantly from their exact values in the case of few-nucleon 
systems. Further reducing the values of the cutoff A leads 
to a significant suppression of the repulsive part of the 
interaction which becomes visible. In another words, if 
one would be able to extend the unitary transformation 
formalism to more-than-two-nucleon Hilbert space, ad- 
ditional three- and morc-nucleon forces would be gen- 
erated through eliminating the high-momentum compo- 
nents. These forces would probably be repulsive [2lJ (at 
least for the systems discussed in this work) and would 
restore the values of the few- and many-nucleon binding 
energies and other observables found in calculations with 
the original, not transformed NN forces. 

On the other hand, it is well known that all realistic 
NN interactions underbind light nuclei such as 3 H and 
4 He. As a consequence, attractive three-nucleon forces 
are needed in order to reproduce the experimental num- 
bers. It is, therefore, in principle possible, that the re- 
pulsive effective many-nucleon forces generated through 
elimination of the high-momentum components compen- 
sate, to some extent, these missing attractive forces, min- 
imizing the total effect of the many-nucleon interactions. 
If this oversimplified picture reflects the real situation, 
one would observe a better agreement with the data for 
calculations based on the LMNN interactions as com- 
pared to the ones based on realistic NN forces. Although 
the results provided by existing many-body calculations 
as well as by the present few-body calculations with 
LMNN interactions appear to be quite promising in this 
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respect |22| , and thus might speak in favor of the above- 
mentioned tendency of counteracting these two kinds of 
many-body forces, to the best of our knowledge, no gen- 
eral proof of the validity of the above-mentioned picture 
has yet been offered. 

In addition, one should keep in mind that only a re- 
stricted information about many-nucleon forces, which 
might have very complicated spin and spatial structure, 
is provided by the discrete spectrum. A much better 
testing ground is served by a variety of scattering ob- 
servables. It would therefore be interesting to test the 
LMNN forces in the three- nucleon continuum. Further, 
we notice that the choice of the LMNN is not unique. 
Preserving the half-on-shell T matrix does not seem to 
be advantageous in any respect, at least from the con- 
ceptual point of view. If this condition is not required, 
infinitely many equivalent LMNN can be constructed by 
means of a unitary transformation in the low-momentum 
subspace. Even if many-body effects appear to be small 
for one particular choice of LMNN forces, this might not 
be the case for a different choice. 

Certainly, the LMNN interactions with small values 
of the cutoff, A ~ 2 fm -1 , are not disadvantageous com- 
pared to existing realistic NN forces regarding the nuclear 
structure calculations, if these realistic forces are viewed 
as purely phcnomcnological parametrizations with no 
physical content, underlying the only requirement of re- 
producing properly the low-energy NN data. 

IV. SUMMARY 

The LMNN interactions have been derived through 
a unitary-transformation theory from realistic nucleon- 
nucleon interactions such as the CD-Bonn and the Nijm-I 
potentials. We have constructed the LMNN interactions 
by two different methods which are based on the com- 
mon unitary transformation. In order to have a cross 
check of both computer codes, we have shown that the 
LMNN interactions obtained by the two methods yield 
the same results. The LMNN interaction reproduces the 
low-energy observables in the two-nucleon system with 
high precision, which has been confirmed by the calcula- 



tion of the deuteron binding energy and the phase shifts. 

The LMNN interaction has been successfully applied 
to the Faddeev-Yakubovsky calculations for three- and 
four- nucleon systems. The calculated binding energies of 
the few-nucleon systems begin to deviate from the values 
calculated using the original NN potentials for A smaller 
than ~ 5 fm -1 , whereas the results obtained by simply 
cutting off the high momentum components without per- 
forming a unitary transformation deviate considerably 
even at much higher values of A. In an appropriately 
truncated (i.e. with A > 5 fm -1 ) low-momentum space 
the LMNN interaction reproduces the exact values of the 
binding energies, at least, for the few-nucleon systems. 

However, we should keep in mind that the calculations 
of ground-state energies using the LMNN interaction for 
A ~ 2 fm -1 yield considerably more attractive results 
than the exact values. We note that as shown in Fig. 1 
of Kuckei et a/.'s work Q one needs more than 4.0 fm -1 
as the cutoff value A in order to reproduce at least quali- 
tatively the saturation property of nuclear matter. Thus, 
the application of the LMNN interaction to structure cal- 
culations should be done with care, though the LMNN 
interaction for A ~ 2 fm -1 may be suitable for the calcu- 
lation of the excitation spectra of low-lying states as has 
been shown in the shell- model calculations H ■ 
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